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J> ■ Abstract. The structure of a light cone in the Godel universe is studied. We derive 

. the intrinsic cone metric, calculate the rotation coefficients of the ray congruence 

forming the cone, determine local differential invariants up to second order, describe the 
crossover (keel) singularities and give a first discussion of its focal points. Contrary 
to many rotation coefficients, some inner differential invariants attain simple finite 
standard values at focal singularities. 
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1. Introduction 

Godel's rotating cosmological model [12], [13] is one of the most interesting solutions 
of Einstein's field equations with negative A-constant, particularly in view of its 
contribution to our understanding of rotation in relativity and its signs of causality 
breakdown due to the existence of closed timelike curves [IB], [2Zj, |2T] - 

The Godel solution is also of interest for a study of light ray caustics, which are basic 
for a discussion of the strong lensing effects in the Universe [2E] , [2H] • There are now many 
papers discussing singularities on characteristic manifolds of the Einstein field equations, 
mainly based on powerful mathematical theorems of Lagrangian and Legendrian maps 
[%] . [TH] . |3] . [H], [II]- One may also mention an older article by Laurent, Rosquist and 
Sviestins jTS] . where the cone of an Ozsvath class III metric [2IIJ was studied. This 
metric already includes the Godel metric as a particular case. 

Focal subsets (caustics) are likely to be present on null hypersurfaces, if a weak 
energy condition holds for its lightlike generators (see [22], [0]), hence almost always 
in realistic astrophysical or cosmological situations. They have often a complicated 
structure and are an obstacle for attempts to solve the characteristic initial value 
problem for Einstein's field equations numerically, since integration along null geodesies 
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runs into difficulties at caustics [Sj, |HI|,[2H]- It would be extremely helpful if existent 
algorithms could be modified or replaced to allow numerical processing through such 
singularities. A preliminary step is to study caustics in exact solutions of the field 
equations. 

In this connection the Godel light cone could be a useful object, since here the 
behaviour of light rays is already sufficiently complex to give an impression of features 
which we can expect in more realistic geometries, and on the other side it is simple 
enough to allow a complete analytical treatment. Due to the five-dimensional group 
of isometries admitted by the Godel metric, which has a four-dimensional transitive 
subgroup, all light cones have the same internal structure. 

A first discussion of the inner geometry of the Godel cone was given by us in 
1972 P, j2], based on the integration of geodesies performed by Kundt in 1956 [T7j . 
Use of computer based formula manipulation technique has shown, that a number of 
complicated relations can be simplified considerably. In particular, the structure of 
caustics and the resulting startling cyclic lens effects found in [T] and j2] now became 
more transparent. The lens effects arise from a quasi-periodic re- focussing of the 
generators and are surprisingly similar to those discussed by Ozsvath and Schiicking 
for the light cone of a plane gravitational wave propagating in vacuum, one of their 
anti-Mach metrics |19j . 

We apply the geometrical decription of null hypersurfaces developped in [Sj,[Z!- 
After integrating the null geodesies in Section II, we derive the intrinsic metric of the 
Godel cone in Section III, calculate and discuss its rotation coefficients and differential 
invariants in Section IV and turn to a description of caustics in Section V. While rotation 
coefficients of the light ray congruence forming the cone have as a rule singularities on 
focal surfaces or keel points, some local inner differential invariants have simple finite 
limits there. It is interesting that this feature - with the same asymptotic values of 
invariants at singularitites - has shown up in all nontrivial light cones studied so far by 
us. A method to relate intrinsic cone coordinates to the angles (6, 0) on the observer 
sky is described in an Appendix. 

2. Light rays in the Godel universe 

2.1. General Congruence 

Godel's stationary solution of Einstein's field equations with cosmological constant 
describes the gravitational field of a uniform distribution of rotating dust matter, where 
- loosely speaking - the gravitational attraction of matter and the added attractice force 
of a negative A-constant is compensated by the centrifugal force of rotation. Hawking 
and Ellis [T5] introduce the Godel metric with the coordinates t, x, y, z as (we have 
exchanged x, y to reach conformity with our notation, the signature will be taken as 
(-1,1,1,1), and the conventions of the Misner-Thorne- Wheeler book will be adopted): 

ds 2 = -dt 2 + dy 2 - X -e 2 yl b dx 2 + dz 2 - 2 e v /b dtdx. 
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The matter density is given by Kfi = -2A = 1/b 2 . With y = y/2be~y/ b one obtains the 
form of the metric used in [2] and employed also here: 

ds 2 = -{dt + —dxf + -(dx 2 + dy 2 ) + dz 2 (1) 

y y 2 

As first shown by Kundt [T7] . the differential equations x" M + x^x"' = for the 
geodesies admit the first integrals 

t > = x (°)> =(-c 2 /V2 + V2y)/c , (2) 

x' = xU' = y(c 2 -y)/(bco), (3) 

y' = xW' = y(x- Cl )/(bc ), (4) 

z' = x (3)/ = c 3 /c , (5) 

Here the prime denotes the derivative with respect to a running parameter s on the 
geodesic, i.e., s is the proper time or invariant length for a non-null geodesic and an 
affine parameter for light rays. The integrals depend on four parameters Cq, c\, c 2 , C3 and 
are subject to the normalization condition g liv x iil x vl = constant, where the constant 
is -1 for a timelike, 1 for a spacelike and for a null geodesic. For the null geodesies 
discussed here, different values of cq correspond only to different definitions of the affine 
parameter, so we assume Cq = 1 subsequently. With the normalization condition 

becomes for null geodesies 

(x - Cl f + (y - c 2 ) 2 = \c\~ 4 = cl (6) 

thus the projection of null geodesies into the xy-pl&ne are confined to a circle with radius 
I I . Light can move to arbitrary large distances only in the z-direction, the direction 
of the rotation axis. Equation (JHJ) is solved with 

x = ci + C4 sin y = c 2 + C4 cos $, (7) 

where $(s) is an unknown function. Since y > 0, no circle point can lie below the y-axis. 
The parameter range for null geodesies is thus constrained by 

c 2 > \/2> 4 |, (8) 

corresponding to points below the line ACD in the upper half plane and above the line 
AEF in the lower half plane of Fig. 1. 

© and (HJ) lead to a single differential equation for $ : 

6$' + c 4 cos $ + c 2 = 0. (9) 

We first shortly discuss the particular case C4 = 0, when the circle in the xy-plane 
shrinks to a point. Here ©-© have the solution 

c 2 c 2 

t = -j=s + t , x= ci, y = c 2 , z = -j=s + z , (10) 

thus through every point P of the Godel universe passes one exceptional light ray: It is 
the light ray sent or received by a comoving observer at P in or opposite to the direction 
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of the local rotation axis. Returning to the general case, integration of JSJ) gives with 
an integration constant k\ 



tan- = -W tan & 2 - c|— — . 11 

2 y c 2 — c 4 v 2o 

The inequality (jHJ) ensures that the roots are real. To simplify the representation, we 
define a new real constant k 



* = ^v / 5 r d, (12) 

and introduce a new affine parameter u> instead of s: 



s 



w = 7^\/4-ci (13) 

Integrating also the remaining equations J2J,©, one finally obtains (e = ±1, since C3 
can have both signs) 

t(w) = 2y / 26arctan (J^tEi tan (w + k)) - + c ^ (14) 



C 2 - C 4 - C j 

cJc 2 2 -cl sin(2w + 2fc) 

x(w) = c\ — — , (15) 

C2 — C4 cos (2w + 2k) 

r 2 - r 2 

yf W ) = ^— ^ -, (16) 

c 2 — c 4 cos (2w + 2k) 
\/c 2 - 2c 2 

z(w) = V2ebw V * = + c 6 (17) 

\jc\-c\ 

as parameter representation of the null geodesies. Counting the number of independent 
parameter one sees that (fT4^) - (fT7|) is the generic null congruence of the Godel cosmos. 
Its explicit form helps to answer questions on null geodesies in the Godel cosmos. For 
example, one can easily conclude that there are no closed null geodesies, which would 
require x^{w) = x fJ, (wi) for some values w and W\. Taking first /x = 3, Equation (|T7jl 
shows that c 3 = or c 2 = v2c4 is needed, which corresponds to z = const. An inspection 
of the relation for \i = 1,2 shows, that these relations can be satisfied by means of 
periodic functions. Thus a subset of null geodesies may return to the same space point, 
but the point is (repeatedly) reached at different times t, due to the aperiodic term 
proportional to w in (JHJ). 



2.2. Light cone geodesies 

We are here interested in those null geodesies, which form a cone with the vertex at 
a point P with the coordinates (0,0,6,0), say. P corresponds to the origin of the 
Hawking-Ellis coordinates. Since all light cones of the Godel cosmos have the same 
intrinsic structure, we could have chosen any other origin in principle. Furthermore, for 
definiteness, the past cone (w > 0) will be considered. Assuming w = at the vertex, 
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we have four relations which will be used to determine ci, C5, k and Cq in terms of the 
remaining parameters C2 and C4: 



Cl 



c 4\/ c 2 — c 2 sin2/i; 
c 2 — c 4 cos 2k 



C5 = — 2b^[2 arctan (4 / ° 4 tan k) , (19) 

V c 2 - c 4 

cos2fc= C ^ + 6C2 , (20) 
0C4 

c 6 = 0. (21) 

The requirement that the Ci,c§, k and C6 exist and are real restricts c 2 and C4 beyond 
(JHJ. In PI a pair (it, u) of transversal parameters was introduced to replace (02,04): 

2 c 4 + c 2 - 6 2 c 2 + c 4 

u = -r, v = . 22 

C4 - c 2 + b c 2 - c 4 



Inverting, we have 



6(l + w 2 )(t; 2 + 1) &(l+« 2 )(*; 2 -l) 



° 2 2(w 2 + w 2 ) ' ° 4 2(« 2 + w 2 ) ' 

The map (02,04) — > (w, d) is not everywhere regular, since the functional 
determinant #^ - = -2uvb 2 {u2+l ^ v l~ l) has for real u,v zeros at u = 0,v = 

0, v 2 = 1. The first two arise from using squares on the lhs of (}2*2*|) . the singularity v 2 = 1 
or C4 = corresponds to the exceptional ray introduced above. In terms of u and v, fl20D 
can be written cos 2/l — l? 2_|_ u 2 • It 

is seen from this equation, that u 2 cannot be negative 
for geodesies forming the Po-cone, this also applies to v 2 because of (JHJ). Thus u and v 
as introduced by (|22p are real. 

For more information we refer to Fig. 1. In this parameter plane the points A 
through F correspond to coordinate pairs (c 2 , c 4 ) given by 

^(0,0), B = K1,0), £7 = ^,-^), D^^-A-), 

B = *4±*,- _L_ ), f = k- ^ 1 < 



V2 + 1' v^+l V2-1' V2-1 

Allowed parameters C2, C4 for generators are subject to (JHJ), and lie above the line APF 
as well as below the line ACD. v = const is the equation of straight lines ("parallels") 
starting at A, with v ranging from v — 1 (the C2-axis) to v — 1 + (the line ACD) in 
the upper half plane. For v < 1 the lines w = const lie in the lower half plane C4 < 0, 
ranging from v — 1 through v = 1/(1 + v2) (the line AEF). Curves with u = const are 
straight lines ("meridians") through the point B, ranging from u — (-BC) to w — ► 00 
(BD) in the upper half plane. We shall find it appropriate (Appendix B) to take u 
negative in one hemisphere. As discussed subsequently, the cone generators cover the 
two shadowed triangles in Fig 1, which correspond to certain quadrants on the observer 
sky, e.g., the lines EF and CD form part of the equator. Two copies of the figure are 
required to cover the full observer sphere. For details we refer to Appendix B. 



Godel Light Cone 



6 




Figure 1. 

C2-C4-parameter plane for null geodesies in the Godel universe (6 = 1 assumed). Null geodesies 
have their integration constants C2 and C4 between the lines ACD and AEF. Those forming 
the (past) light cone through Pq are confined to the shadowed triangles, with C4 > for the 
western hemisphere and C4 < for the eastern hemisphere of the observer sky. But only 
the northern hemisphere is mapped 1:1 to the two triangle regions, for the missing southern 
hemisphere a second copy of the figure is needed. To obtain the topology of a sphere, certain 
boundary lines of the triangles in both copies must be identified pointwise. The lines CD and 
EF correspond to parts of the equator, the pole ray at B is the exceptional ray towards the 
rotation direction. 



Returning to the cone representation, the parameters C5, c±, k and c§ can be written 
in terms of u and v in a compact form as solutions of (fT% )) -(f2"I j) : 

u(l-v 2 ) 

Cl = ~ b 2,2 ' 

C5 = — 2v / 26arctanu, 

7 u 
tan k = — , 

v 

c 6 =0. 

Substituting these values into (0)-©, one obtains as a parameter representation of the 
light cone through P : 

1 \ bw 1 r- (u 2 + f 2 )tanw 

t = t(u,v,w) = -={v + -) + 2V26arctan — — — — — , (23) 

y/2 v v{u 2 + 1) + u(v 2 - 1) tan(w) 

, ,, 9 . sin w(v (u 2 — 1) cosw + u(v 2 + 1) sinw) 
x = x{u, v, w) = b(v 2 - 1 ^ ■ \ 2 , 2l ' ■ \ 2 , (24) 

f s bv 2 (u 2 + l) 

y = y{u, v, w) = : ■ : 25 

(v cos w — u sin wy + v A (u cos w + v sin wj z 

6eif 



z = z(u, v, w) = — j ^\/ < ov 2 — 1 — v 4 . (26) 
V2f 
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Figure 2. 

The twisted generators of the past light cone through Pq are shown in a fictitious Euklidean 
3-space 1Z 3 . Plotted are the coordinates x(u,v,w),y(u,v,w),t(u,v,w) from (|23 |) -l|25 |) for rays 
on the parallel v = 1.5 and several u with constant separation 7r/30 in the meridian angle 
4> = 2arctanu corresponding to u. The z-coordinate is suppressed, keels appear as points, 
focal surfaces as lines. The affine parameter ranges from w = at the vertex (top) tow = 1.57T 
(bottom), the plot ends before the second keel line at w = 2tt is reached.- Distances cannot 
be represented correctly in such a projection, the type of caustics is preserved however, since 
the map from the curved V s to 7Z? is a diffeomorphism. 



Positive values of the affine parameter w correspond to the past light cone, negative 
values to the future cone. The sign e distinguishes between the northern (e = 1) and 
southern (e = —1) hemisphere of the observer sky. We note the following invariance 
property of the system (J2S1)-(|2E1): Substituting — 1/u for u and 1/v for v leads to the 
same geodesies: 

x^(u,v,w) = x t *(--,-,w) (27) 

U V 

The same map sends also (c 2 ,c 4 ) into (c 2 , — c 4 ). Some geodesies of the parallel v = 1.5 
are plotted in Fig. 2. 
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Knowing the tangential vector ^p- of the past light cone from (J23j) - (j26|) or from 
|-©) one can calculate the redshift 2; of distant objects from the well-known relation 

-, , V " ijavj emitter 

1 + z - 



\ dm V V Q nv) observer 

With V^ 1 = 5q and (|23jl - <|2fijl one finds z = 0: As noted already by Godel in his original 
paper jT^j, distant objects comoving with the cosmic fluid would show no redshift, 
proving that this model cannot represent the real Universe. 

3. Light cone metric 

The spacetime coordinates of light-rays through P , (f2~3"j) - (f2T)J) . depend on the afline 
parameter w as well as on two quantities u, v. While w determines a position on a light 
ray, u and v label a ray. The tripel w l = (w,u,v) may therefore be used as intrinsic 
coordinates on the cone. Below we will see how u and v are related to the angles 8, <f) 
on the sky of a comoving observer at Pq, who wants to fix an event on his past light 
cone. The intrinsic three-dimensional metric of the light cone at Pq can be found from 
(pS Jl -pfi jl by means of 

dx^ dx 11 

lik 9uu~n 7T{ Zi 

where w % = (w,v,u). Since ^ is the tangential vector to the cone and hence null, the 
components of the cone metric can be reduced to a two-dimensional metric ^ab'- 

7oo = 0, (28) 
7oa = 0, (29) 
b 2 

722 = ~TT ( w2 f? + 2w hfi sin w cos w + fl sin2 w + hfl sin 4 w), (30) 

b 2 ff sin 2 w(w — sin w cos w) 

723 = 4 , 2 7T , (31) 

f 4 (M^ + 1) 



t/ J X Olll U/ V -±U — J X Dill LUJ 

^ = ^J77a~ZT\2 > ( 32 ) 



b 2 fl sin 2 w(Av 2 - f 2 sin 2 w) 
v\u 2 + If 

(33) 

where we have introduced the two functions 

f 1 (v)=v 2 -l, f 2 ( v )=Qv 2 -l- v \ (34) 

which are useful to compactify expressions. 

In general, the ^ab form an one-dimensional sequence of positive-definite two- 
dimensional metrics on the cone, parametrized by the afline parameter w. They 
represent metric spheres in the neighbourhood of the vertex w — 0, but become 
progressively deformed for increasing parameters w. On some subsets of the cone the 
inner metric degenerates additionally, i.e. the two-dimensional determinant 7 becomes 
zero. This signifies an intersection of light rays forming the cone, if the zero of 7 
is not caused by a coordinate singularity. In focal points or caustics geodesies with 
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infinitesimally differing values of u, v meet, while in crossover or keel points (we have 
taken the latter notation from a paper by M. Riesz |26j ) the intersecting geodesies may 
have quite different transversal parameters. Sets of focal points form in general two- 
dimensional focal surfaces on a null hypersurface. For the Godel cone the determinant 
7 can be represented in a compact form: 

7 = \iab\ = h 2 , (35) 
2b 2 fip(v,w) sin w 

(U 2 + l)V/2^ 

p(v, w) = f 2 w cos w + f 2 sin w (37) 

(abbreviations are deliberately chosen in this paper to compactify expressions). The 
simple expression for 7 allows to pick up the cone singularities easily. Apart from the 
vertex w = 0, 7 vanishes periodically for w = rnr, n an integer, which is similar to the 
behaviour of light cones in a closed Robert son- Walker model with the time extended 
to several cycles. But contrary to the RW case, the two-dimensional spacelike surfaces 
w = const do not shrink to points at w = rnr, n =/ 0, but rather to spacelike lines, the 
keel lines, discussed in Section V. Further zeros of 7 are given by p(v,w) = 0, which is 
the equation of the focal surfaces, also discussed in Section V. The remaining zeros of 7 
are given by fx — or v — 1, corresponding to the exceptional pole rays, and u —>■ 00. 
The latter is related to a coordinate singularity. 

For later use it is appropriate to introduce a second function q(v,w), with the 
property of being not negative: 

q(v,w) =4(1 + fx) -f 2 sin 2 w. (38) 

q vanishes if and only if f'2/f 2 — —cos 2 w. Since the lhs is not negative, the condition 
q = holds only for cosw = (or w — (1 + 2n)7r/2, n integer), and for f 2 = 0, or 
v — 1 + y/2. Since at these points also p = 0, the equation q = represents curves on 
the light cone, where the equator rays v = 1 + v^2 meet the nth focal surface (Section 
IV). With q, the intrinsic cone metric can also be written as 

b 2 

722 = -r-, ^— (P 2 - 2p/ 2 sin 3 w + qf 2 sin 4 w), (39) 

v J 2 cos w 

b 2 fx sin 2 w 

723 = 4/ 2 1 -n (p-qsmw), (40) 

iriir + 1) cos w 



^ 33 = ^,2_ L n 2 g» ( 41 ) 



b 2 f 2 sin 2 w 
v 4 (u 2 + 1) 

allowing to check easily that the rank of 7^ indeed becomes 1 (i.e. |7ab| = 0, but not 
all 'Jab = 0) at caustics p = 0. 

The range and the geometrical meaning of the transversal cone coordinates u and 
v must now be discussed. Apparently, u may take all values out of the range (— 00, 00). 
For a real z-coordinate the function f 2 = Qv 2 — 1 — v 4 cannot have negative values, so v 
is restricted to an interval (v min , v max ) where v m i n = 1/(1 + y/2) and v max = 1 + y/2. If 
v belongs to this range, also 1/v does. However, the substitution (|2*7j) shows that two 
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different pairs (u, v) might represent the same null geodesic. Also two different geodesies 
might be represented by the same (u, v )-pair. 

A natural way to parametrize a light cone is to take the angular coordinates (9, (f) 
at the sky of an observer sitting at the vertex and comoving with the fluid. As usual, 9 
ranges from (north pole) to 7r (south pole), and <ft from to 2tc. To find the relation 
between the sky coordinates and (u, v) we have used a method, which is exclusively 
based on the intrinsic cone metric. Leaving the details for Appendix B, the result is 

, 1— cos 6 9 v^ + sin^ , 1 — u 2 . . r-v 2 — 1 , ,„ s 

u =TT 1> v =X 7n ~ai cos0=-— -, sm9 = V2^—. 42 

1 + cos <p V2 - sin 9 1 -f- u v l + 1 

u ranges from u = for <p = to oo at <fi = ix, jumps there to — oo, and increases to zero 

at = 2ix. v starts from v = 1 at the north pole (9 = 0), increases to v max = 1 + v2 

at the equator (9 = 7r/2), and decreases to 1 at the south pole (9 = n). As explained 

in the Appendix, only the partial interval (l,v max ) is used for v. A point on the sphere 

is fixed by a pair (u, v) together with the sign of e. This ensures that a pair (u, v) from 

the w-range (— oo, oo) and t>-range (1, v max ) is a one-two-one map of the light rays in the 

northern resp. southern hemisphere. The exceptional ray with v = l,e = 1 corresponds 

to the north pole, the antipodal ray with v = 1, e = — 1 to the south pole. 

Transforming the inner metric 7^ to angular coordinates 9, <fi does not simplify 
neither its form nor other relations very much, so we continue to work with u and v as 
transversal coordinates. 

In the chosen representation the equator is given by v = v max = 1 + v2 and 
geodesies sent out in these directions (orthogonal to the rotation axis) lie always in the 
plane z = 0. Since here ji = 0, the cone metric becomes singular, its determinant 7 
tends to infinity. This is a coordinate singularity: The angular coordinates are regular 
along the equator, but the functional determinant |§^a| suffers from a diverging factor 

— 1/2 

/ 2 .If carefully treated, this divergence will not cause trouble. 



4. Rotation coefficients and invariants 

4-1. Geometries on null hypersurfaces 

The local differential geometry of null hypersurfaces such as a cone was in some detail 
described in and [7j, see also [211,123 • We summarize the points most important for 
us. This geometry is formulated in terms of the rotation coefficients of a certain class 
of triads, defined as follows: At every regular point P of the cone there exists a unique 
direction e\ i = 1,2, 3, the direction of the null geodesic passing that point. e l satisfies 
7^e fc = and is given up to a factor by e l = S\ in the coordinate system (u,v,w) used 
in the last section. The two other directions, which are spacelike in regular points and 
orthogonal to each other, may be combined linearly to form a complex vector t l . We 
have t t e k, -fik = and normalize f such that t l t kr )ik = 0, t*t fc 7ifc = 1. The transversal 
directions t l are determined only up to a transformation t n = e lu (t % — Re 1 ), u real and 
k complex. Note that e l is also subject to a change e' 1 = \e l (A real), since the running 
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parameter along a ray may be chosen arbitrarily (it need not to be an affine parameter). 
The covariant components of the transversal directions are given by tj = 7^^ and ji, 
where 7$ is defined by t % ^i = 0, e*7i = 1. This completes the covariant triad. The 
rotation coefficients divergence p, shear a as well as other coefficients are given in terms 
of the derivatives of the triad: 



p + IV 
a 

T 

X 



£ l i k (ii,k — tk,i), 

t t (ti t k l-k.i) ; 

^?e fc (7*,fc-7M), 
^ fe (7i,fc-7fc,i)- 



(43) 
(44) 
(45) 

(46) 
(47) 

t is related to the intrinsic geometry of the two-dimensional wave surfaces w = const, 
with w here as an affine parameter of the generating null geodesies. The coefficients x 
and if reflect properties of the triad, which are geometrically not relevant. In particular, 
if 7, is chosen as gradient, as we will do here for simplicity, both x an d <f are zero. A 
change of the triad 



1 

— 7i ~\~ K>H + fttj, 



e tuJ (f - R\e l 



produces a change of the rotation coefficients as follows: 



p 


= Xp, 






a' 


= \e- 2iw o, 






T 




inXDuj - 


inXv + k\o~ — nXp), 


v' 


= X(v + Du), 






x' 


= 5«^ + T 


— np\ + 


ivnX — RXo~ + XDk), 




1 ■ !/ 


Kf + 5k - 




iip' 


= -t<p+-{KT- 


-8k) 



+ 



sx 
T 



inXv — kXct + XDk) 



k(2x + -r- + zkXv — kXo~ + XDk). 

2 A 

The last two equations show that x = an d if = are preserved for k, X satisfying 
5k — 5k + kt — kt = 0, 5X/X 2 + Dk — K(p + iv) — kg = 0. Coordinate invariant 
statements are formulated in terms of those functions of the rotation coefficients and 
their derivatives, which are invariant with respect to the allowed transformations of 
the triad. The group of allowed triad transformations defines the type of null surface 
geometry in the spirit of Felix Klein's "Erlangen program" [TB] . The most important 
geometries are the just outlined inner geometry and the affine geometry, where the 
concept of an affine parameter for the rays is given as additional geometrical element. 
For other geometries on null hypersurfaces and for more details we refer to or [7j, 



Godel Light Cone 



12 



and for a similar definition of null surface geometries the papers by Penrose [21], [2H] 
should be consulted. 

4-2. Application to the Godel cone 

We first determine the rotation coefficients for the Godel cone. The divergence is 
calculated from p = — and may be written as 

p = -cot2w , (48) 

2p cos w 

with the functions p(v,w),q(v,w) defined by (jH7j) . (jHHj) . The divergence tends to ±oo 
at w = nit (keels) and p = (focal surfaces) and becomes zero at the two-dimensional 
surfaces 

tan2w f? , , 

-— ~i m 

between focal surfaces and keels: Since the rhs of this equation is not negative, the 
range of w, where p = is possible, and hence the position of a zero-divergence surface, 
is restricted by the condition 

(2m- 1)- < w < m-, m= 1,2,3.... (50) 



The amount of shear follows most easily from another general relation \o~\ = p — 

dw 



2 q 2 8v 2 cos 2 w + q(2 sin 2 w - 3) 1 ff 

°1 — ~r^> o 1 ^ — : 5 !~ A ■ 2 5 ~r^- 

Ap z cos z w 2p sin w cos z w 4 sin w cos z w Av z 

Like the divergence, also the shear goes to infinity at focal surfaces and keels. The 
remaining nonvanishing coefficients may be calculated from ()43|) - (j43)) . the details are 
given in Appendix A. Splitting a into real and imaginary parts, a = a\ + io"2, one 
obtains 

* 1= _^_ + !±^ + g ; (52) 

sin 2w 2q 2p cos w 
°"2 = , (53) 



2vq 
2vq 

ffciv 2 + l){2v 2 - f 2 sin 2 w) 



(54) 



bf lP VTq ■ (55) 



The triad has been chosen so that the resulting rotation coefficients look as simple 
as possible. There exists in general one (and only one) first-order inner invariant of a 
null hypersurface, i.e. an invariant function formed from the rotation coefficients alone, 
without derivatives. This is the quantity j = A? or any function of j. It is useful to 
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consider 1/j 2 , which measures the anisotropic behaviour of the generators around a 
given one: 

2 pf?(4:V 2 sin w — p) sin 2 w cos 2 w + v 2 (p — q sin w) 2 

= 21 in ■ 2 T\ : V2 • ( 56 ) 

v z (p{2 sm w — l) — qsmw) z 

At caustics p = (and keel points with w = mr) this gives j 2 — 1 or j — ±1, 
depending on the sign of p. Along the two exceptional rays (pole rays) the shear vanishes 
and the anisotropy measure 1/j 2 is zero, in accordance with the symmetry properties 
of the cone. 

From the rotation coefficients and their derivatives one may also form some second- 
order invariants. The purely transversal projections of the four dimensional Ricci and 
Weyl tensor into the cone are closely related to them. The Ricci and Weyl tensor 
projections are -without any addition of further embedding quantity- equal to similar 
projections related to an intrinsic Riemann tensor -R ifc ; m of the null hypersurface, 



MjRuFtm = to = Dp- p 2 -aa, (57) 
tVe^ ™t m = ip = Da- 2a(p - iv). (58) 



A straightforward calculation gives 
(v 2 _ + V 2 
Av 



- = ^5FS (59) 



2 Acos 2 w 1 JJ 2 sinwcosw 
q 2v 2 vq 



To calculate an intrinsic curvature tensor of a null hypersurface as in ([57|58|h one 
needs an affine connexion in spite of the missing unique contravariant metric tensor. We 
again refer to [B],[I| and note only shortly, that in general the resulting affine connexion 
and hence the inner Riemann tensor depend on the triad. Independence holds only 
for special projections such as (|57|58jl . u, ip are nevertheless not yet affine or inner 
invariants of the cone, they are densities, and one has to apply suitable factors of p 
or |er| to generate invariants of the affine geometry. To obtain invariants of the inner 
geometry, one has to take a certain linear combination of these affine invariants. If we 
define 

, uj ib 1 , if Dp Da\ 2v , . 

I = I 1 + tI 2 = i — --^L + --j)= l-H + 61 

p\a\ a\a\ j \a\ \ p a J \a\ 

then this quantity is a second-order differential invariant of the inner geometry. A short 
calculation gives for the real part Ii = (Ds + 2v)j\o\ (where s is the argument of a such 



\a\ 3 h = J1VJ2V L7 (62) 



fl^fhi^v 2 sin w — p) 
Apv ? - 

Ii is a measure of the rotation of the two shear directions (defined as directions 
to neighbour rays with extremal distance change, [U]) with regard to the generator 



Godel Light Cone 



14 



congruence. The imaginary part, which can also be written I 2 = Dj/p, is slightly more 
complicated and may be represented as 

-,2 j_ „• ^3 



, |3j 2 = ?p + lip + l 2 p* + t 3 P , g3 ^ 

Ap 2 v 2 sin 2 w cos w(o sin w + p{\ — 2 sin 2 w)) 



with 



zo = Av 2 q 2 (q — 2v 
i x = ff( v 2 + l) 2 sin 7 w - 8v 2 tf(v 2 - 3)(3v 2 - 1) sin 5 w 

+ 64w 4 (/i - v 2 ) sin 3 w - QAv 6 sin w, (64) 
i 2 = - 2f 2 f 2 sin 4 w + I2v 2 f 2 sin 2 w + 8v A (l - 2 sin 2 w), 
k = - / 2 sinw. 

The Gaussian curvature = — 2rf+5r+5f (see jH]) of the two-dimensional surfaces 
w = const is in general not an invariant of null hypersurfaces. The only exception 
are Killing horizons (sometimes called "totally geodesic null hypersurfaces", see, e.g., 
Hajicek ^3]), defined by the condition that the inner metric admits a Killing symmetry 
with the generators as Killing vectors. In all other cases K depends on the chosen 
foliation and is not significant for the cone geometry: A change of the affine parameter 
as w = a(u, v)w (keeping the vertex at w = 0) leads to a different curvature K. K is 
only invariant under transformations of the transversal parameters x A/ = f A (x B ).Foi 
our foliation w = const an explicit calculation gives 

K=%A (65) 

with 

h = (32v 6 -f? sin 2 w)/(b 2 f 2 ), (66) 

k 2 = Af 2 v 2 {v 2 + l) 2 sinu7(sin 2 W / 2 - 2i; 2 )/(6 2 / 2 ). (67) 

K generally tends to zero for large w, apart from spikes at focal surfaces p = 0. 

In jn], [Zj points on a given ray have been classified according to the focussing 
behaviour of neighbouring geodesies. A point on a ray was called elliptic, if the spatial 
distance to all neighbouring rays either increases or decreases, and hyperbolic, if some 
rays converge and other diverge. The sign of p 2 — \a\ 2 distinguishes both types of points. 
A positive sign (or j 2 > 1 at points with nonvanishing shear \a\) corresponds to elliptic 
points. Evidently, near the vertex at w = all points on all rays are elliptic, this is 
also seen from an expansion of j near the vertex, j w —6v 2 /(f 2 w 2 ). Zeros or infinities 
of p 2 — | a | 2 along a ray may signify the transition from elliptic to hyperbolic points (or 
vice versa, a point with p 2 — \ o~\ 2 = will be called parabolic). It is not difficult to verify 
that p 2 — \a\ 2 can be written as 

p sin w Av 2 

Thus caustic singularities (p = or w = rnc) can be transition points on a ray. We 
illustrate this for the two pole rays {y = 1) and for the equator rays (v = 1 + V2). For 
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pole rays p 2 — \<r\ 2 = cot 2 w, all ray points are elliptic, with the exception of isolated 
parabolic points at w = mr + 7r/2, hence no proper transition point exists. In the case 
of equator rays we have p 2 — \a\ 2 = 2cot2u>/u>, here all points with w = n(2n — l)/4, 
n integer, are transition points, including also non-caustic points. 

To have some idea of the cone structure, we follow a typical ray with 1 < v < 
1 + \/2, u arbitrary, from the vertex w = down the cone (we consider the past 
light cone, where an increasing affine parameter w means decreasing time). Near the 
vertex, the cone resembles the Minkowski light cone with p = —1/w, zero shear and 
exclusively elliptic points. All neighbour rays recede from the chosen ray. Along the 
ray, \p\ decreases from p = — oo at w = and \a\ increases from zero, until a first 
transition point is reached, where both are equal. At the transition point the invariant 
j has increased from — oo at w = to j = — 1. Behind this point a domain of hyperbolic 
points begins, where some neighbour rays start to decrease their distance to the chosen 
ray. The next significant point on the ray is the zero- divergence point, where also j 
reaches zero for the first time. After passing this point, p is positive and increases faster 
than \a\. Thus a second transition point is encountered with j = 1, ending the domain 
of hyperbolic points. Behind the transition point a region of elliptic points begins, all 
rays converge towards our ray, preparing for a meeting at the first focal point. At the 
focal point both p and |<r| tend to infinity, their quotient j jumps from 1 to —1. 

Behind the focal point, p increases from large negative values, passing a third zero- 
divergence point in the interval 3ir/4 < w < it, until the first keel point is reached at 
w = 7r. The whole region between focal surface and keel consists again of hyperbolic 
points. Behind the keel point the cycle starts again, with changed positions of the 
transition and focal points relative to the zero-divergence and keel points. 



5. Focal sets 

5.1. Keels 

The most impressive singularities of the Godel cone are those at w = nir (Fig. 2), the keel 
points. They correspond to "points of the first kind" on the light cone of Ozsvath and 
Schiicking's anti-Mach-metric [19 a . At the keel point all rays with equal v and different 
u meet . An one-dimensional set of connected keel points is denoted as keel line. Every 
integer n gives a keel line with the spacetime coordinates 



nirb(v 2 + 1) nnbey/ Qv 2 — 1 — v 4 . . 

tkeel ' ^keel U, Dkeel 0> %keel 7E ■ \ J 

V2f V2f 
Thus keel lines can be considered as circle segments in the pseudo-Euklidean z — t plane 
with a length L increasing with n (Fig. 3). 

Their end points lie on the exceptional ray and its antipode. For n = the keel 
lines shrink to the vertex P , and for n^O the observer world line crosses the keel lines 
only for equator rays v = v max . It is easily checked that the keel lines are spacelike but 
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Figure 3. 

The keel lines w = nir appear as sequence of isolated sphere segments with increasing length, 
if projected into the pseudo-Euklidean plane z — t. Shown are the keel segments for the future 
and past cone (n = 1...6, b = 1). The £-axis is the projection of the observer world line. The 
endpoints of the keel lines lie on the exceptional ray in the northern hemisphere (z > 0) and 
on its antipodal ray in the southern hemisphere (z < 0). These rays are plotted as straight 
lines. The keel endpoints are also the intersection points of the nth keel line with the nth focal 
surface. Every point on a keel line (which is parametrized by v) is an intersection point of all 
generators with different parameters u (and the same v). In particular, the intersection of the 
keel with the t-axis corresponds to equator rays v = v ma x, the only rays which return to the 
observer worldline. 



not geodesic in the Godel geometry. The components of the tangential vector dx^/dv 



are 



(dx»\ mrb ( 2 e(v A - 1) \ 

and its norm is b 2 n 2 7r 2 ff/(v A f2). The first normal of the keel (the binormal does not 
exist) is the timelike unit vector 

< = /o(\ i> 2 + °> °> eV6n 2 -l-n 4 ), (71) 
V2(n 2 — 1) 

the (first) curvature is found as 

4n 4 

h = nbn{v 2 -iy (72) 
it diverges at the two endpoints v — 1 of the keel lines. Here also the otherwise spacelike 
tangential vector degenerates to zero. The invariant total length of a keel segment, 
ranging over the full observer sphere, is however finite and given by 

r rl+V2 ( v 2 _ 1)2 

L — ds = 2mrb / / ; dv w 2.39628 nnb. (73) 

J Ji v 2 \ / 6v 2 — v 4 — 1 
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The divergence of the integrand at the equator v — 1 + v^2 arises from the coordinate 
singularity there, but the integral converges (we could equally well have integrated from 
v = 1/(1 + -\/2) to v — 1 with the same result). 

It is of interest to study the behaviour of rotation coefficients and in particular 
differential invariants near the singularities p = and w = nn. A power series expansion 
around keel points w = nn, n^O, leads to 

1 2v2 

P = - 7T7 V 72 +o((w-nn)), (74) 

2{w — nn) nnff 



1 2f 

2\w — nn\ nitfi 



\a\ = — j- ZTc2^( w ~ n7r ) + °( w ~ n7r )' (^) 



^ 2 



j = — U(w — n-n) H -^(w — nn)) + o((w — nn) ) (76) 

nnff 

for small w — nn, where U(x) is the step function, U(x) = 1 for x > and —1 for x < 0. 
While the ray divergence as well as the shear amount individually have first order poles 
at focal or keel points, their invariant quotient j remains finite, but jumps from +1 to 
— 1, if w increases. For the second order invariants I\, 1% one obtains near keel points: 

q2 



h= — ^f^(-l) n (w - nn) 3 + o((w - nn)f, (77) 

h = —pji-lTiw - nn) + o{{w - nn) 2 ), (78) 
nnff 

thus the complex invariant I = I x + il 2 vanishes here. 
5.2. Focal surfaces 

As already discussed in previous sections, apart from the keel lines also focal singularities 
exist (we shall not discuss coordinate singularities). In geometrical optics a caustic (set of 
focal points) is the locus where the rays have an envelope and the intensity a singularity. 
Here focal point are similiarly defined as points of intersection of infinitesimally close 
geodesies, which satisfy the relations x^{u + 5u, v + 6~v,w + 5w) — x^{u, v, w) = 0, where 
x^(u,v,w) is the cone congruence Expanding, we have 

<9x M , dx^ , dx 11 , , , 

— Su + — 5v + —5w = 79 
ou ov ow 

for suitable displacements 5u, 5v, 8w. Let us assume v ^ 1 (we let out the exceptional 
rays) and sin if ^ (no keel points are considered). We can then eliminate the 
displacements from (|T9*|) and obtain two relations, which must be satisfied for the 
coordinates u, v, w of focal points on the cone: 

p(v, w)r(u, v, w) = 0, (80) 
p(v, w)s(u, v, w) = 0. (81) 

Here r = r i u% an d s = J2 s i u% are polynomials of fourth order in u, r% and Si are 
polynomials in v, w, sin w and cos w. A closer inspection shows that with the restrictions 
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sinw ^ 0,t) / 1, r and s cannot vanish simultaneously, thus we conclude that focal 
points are given by p(v, w) = or 

(v 2 — l) 2 w cos w + (6v 2 — 1 — t> 4 ) sin w = 0, 
which might also be written as 

_ta™ = £ 1 an2 , 

Alternatively, focal points can be considered as the critical points of the map (w, v, w) — > 
a; M , i.e, points where the rank of the Jacobian matrix is not maximal [3]. This leads to 
the same condition (J%2*|) . (jH2J) is the equation of two-dimensional surfaces on the cone. 
It is similar to the condition for "points of the second kind" on the Ozsvath-Schucking 
cone jinj. Note also the similarity of this equation to the equation for zero- divergence 
surfaces f!49j) . As in this case we conclude that focal surfaces can only occur in regions, 
where the affine parameter w is confined to the intervals 

(2n - 1)- < w < mr, n = 1,2.3..., (83) 

since only here — tanw/w is not negative. Thus the focal set decays into an infinite 
number of separated two-dimensional sheets. The circular functions generate a quasi- 
periodic behaviour with similar but not identical shapes for the sheets. 
If we solve ()82)1 for v, we find 

2 2(sinu> + (— l) n+1 \ / 2 sin 2 w — wsinwcosw) 

v — 1 H . (84) 

smw — w cos w 

If w is in the interval (|83|) . the square root is real. Since v 2 — 1 varies only between 
and 2 + 2^/2 w 4.8284, the quotient in (|H4^1 must fit into this interval. This is achieved 
by choosing the sign of the square root as indicated. 

Keels and focal surfaces are not completely separated. Every keel line w = mr has 
two common points with the nth sheet of focal surfaces, corresponding to the keel line 
endpoints w = nn,v = 1. This common point lies on the exceptional null ray resp. its 
ant ip ode, see also Figure 3. 

If (JH1J) is inserted into the equations (f2~3*j) - (f2T)|) . one obtains an explicit representation 
x^ ocal (u, w) = x^{u, v[w], w) of the focal surfaces. The tangential directions at the point 
iu,w) on a focal surface are spanned by the spacelike vector dx^ ocal /du and the vector 

dx^ dv* dx^ dv* dx^ (u 2 + 1) , s 

k = j>- + i~^r + i-^Tf ' 85 

ow aw ov aw ou ji cos w 



where v*(w) is the function defined by (J84J). fc M is in general spacelike on the cone, 

but becomes a null vector at focal surfaces. We expect that the curves to which is 

tangent are non-geodesic null curves (see (2] for an excellent discussion of non-geodesic 

null lines in the Minkowski spacetime). 

Considering the invariants near focal surfaces, it turns out that j jumps from +1 

to —1 with increasing w - this is the same behaviour as in keel points. Near the nth 

focal surface p = we may write with the step function U(p): 

NriTT/ . . 4 cos 2 w( f? sin 2 w — 2v 2 ) . . 9 . . , 

J = (-l) n U(p)(l ^ J -p) + o(p 2 ). 86 

q 2 sin w 
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Note that with increasing w the function p(v, w) reaches zero at the nth focal surface 
from positive (negative) values, if n is odd (even). A similar calculation for the invariants 
h,I 2 gives 

h = 8( ~ 1)W+1/l2v ^ C ° s3wSinW P 2 + o{p\ (87) 
vq 6 

8 cos 2 w(2v 2 — f? sin 2 w) . , 9 . . . 

h = 2 • -b + o(p 2 ) 88 

q z sin w 

near the focal surfaces p = 0. 

This shows that at both focal and keel points the complex invariant I — Ii + il 2 
vanishes, while j is 1, modulo a sign. The same result holds for the Oszvath-Schucking 
light cone [TH], see 0. 
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Appendix A. A triad on the Godel cone 

For an explicit calculation of the rotation coefficients ()43|). ()44jl. (j4"o"|) one needs the 
components of a suitable triad e k ,t k , t k and 7*, tfc, i k on the light cone. We have already 
chosen e l = 8\. The inner metric is given in terms of the triad by 

lik = Uh + Ut k (A.l) 

Comparing this expression with (|25 j) - (|32j) shows that U = 5ft a, t % = S A t A , A = 2, 3. It 
is not difficult to verify that 



bsinwip — qsmw) bp I 2 

*2 = - 2 ts= — + \—r, A -2 

v z cos W\J Lq v v qfe 

= ^qbfrsmw 

3 y/2v 2 {u 2 + l) 1 ' ^ 

reproduces the equations (jHSJ-ljUJ)- For the contravariant components we use the 
normalization conditions t l t k ^fik = 0, t l t kr )ik = 1 and obtain 

t 3 = v2 ( u2 + jj + A u2 + l)v72(gsinw - p) , A g , 
6/i sin w 26/i cos Wy/2qp 
The nonvanishing rotation coefficients may then be found from 

p + iu= - t 2 t 2 ,i - t 3 t 3 ,i, (A.6) 

a = - Pt 2 ,i - Pi 3 ,i, (A.7) 

= (t 2 t 3 -t 3 t 2 )(t2,3-t3 >2 ). (A.8) 
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Appendix B. Observer sky and cone parametrization 

In sky coordinates 9, (f> any cone metric can be expanded in powers of an afflne parameter 
w* near the vertex (see, e.g., 0): 

w* 2 

lee = — + o{w ), 

= o(w* 5 ), 
w* 2 

= rL-wn 2 e + o(w* 4 '). 

A similar expansion of the Godel cone metric in powers of w gives 

^ 2 b 2 f 2 
^ = v 2 (u 2 + l) 2+0{w) > 

723 = 0(W 5 ), 

iQw 2 b 2 4 

733 = — f V o{w ). 

h 

The coordinates (u,v,w) are related to (9,(p,w*), and this coordinate transformation 
should take approximately the form 9 = 9(u,v), <fi = <p(u,v), w* = w/m(u,v) near the 
vertex, i.e., for small w. For the transformation functions we thus obtain the differential 
equations 

fao\ 2 1 (d<p\ 2 zin2e _ mw 

\8u) \8u J v 2 {v? + l) 

aeae a<t>d<j> . 2 

7T7T + sm 9 = °' R2 

8u 8v 8u 8v 

which can easily be solved, if we assume 9 = 9{v),(f) = (f)(u). Then (|B.2|) is already 
satisfied and the other two give 

(w 2 + 1 



<9ti fsin^ <9t> v/2 

where e\ = ±1. The second equation here shows that m depends on v only, the first 
equation then says that both sides must be equal to a constant k\ independent of u and 
v. Integrating, we first obtain 

u = tan (<f)/ki). 

The differential equation for 9 follows as 
1 89 2e x k lV 



sin9 8v j\\Th 
and is solved by 



(B.4) 



n / £ \ eifel/2 
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where k 2 is a second constant. To obtain a real square root of f 2 , v had to be confined 
to the interval v min = 1/(1 + \/2) through v max = 1 + \[2. We refer to (t> m in, 1) as the 
min interval and to (l,v max ) as the max interval. 

To determine the coefficients hi, k 2 in Eq. ()B.5|) . we must face the possibility that 
they differ in different parts of the sphere. 

The symmetry properties suggest that the poles of the observer sphere are related 
to the local rotation axis and are thus given by v — 1. We consider (jB.5|) near the north 
pole and assume, that the v belong to the max range. Expanding the rhs in powers of 
small v — 1, on obtains k 2 ((v — 1) / '2) eifcl//2 , thus the sign of e\k\ must be positive to ensure 
that the rhs vanishes for v — > 1, as does the lhs. We also conclude that k 2 > 0. A similar 
conclusion is reached for the min range of v: Apart from the positive sign of eiki, also 
k 2 (— l) £lfcl / 2 must be positive and real. (Note, the mm-interval of v is obtained from the 
max-interval by applying the map v — > 1/v, the expression within the bracket in (|B.5j) 
attains the factor -1 under this map). Using the identity tan (|) tan (| — |) = 1, one 
can easily repeat the calculation near the south pole. It is seen that the lhs of (|B.5|) 
diverges there, thus also the rhs diverges, and this requires eiki < 0, holding again for 
min as well as for max ranges of v. Furthermore, we have for the max (min) range 
k 2 > (A; 2 (-l)- fc 3/2 > o). 

Moving now from the north pole towards the equator, assuming the max interval, 
v as well as 9 increase until v = v max is reached, which corresponds to a 9 max = 
2 arctan [k 2 /2 eikl ^ 2 ]. 6 max cannot represent the other pole 9 — ir, since the lhs of (jB.5|) 
diverges at 9 = it, while the rhs here is regular. Thus only part of the sky is covered 
by v values in the max range. It is convenient to assume that this part is the northern 
hemisphere, i.e., 9 max = ir/2. This fixes k 2 by k 2 = 2 eifcl / 4 . If we had started our walk 
in the min region of v, taking v — 1 at the north pole and decreasing v to the equator, 
we would have obtained a similar conclusion: Also the min range covers the sphere only 
from the pole to the equator, here k 2 is fixed by k 2 = 2 eifel / 4 (— l) _eifcl / 2 . Our walk could 
have started from the south pole, reaching the equator from the south, but the results 
for k 2 are the same. 

Further conclusions depend on k\. Mapping the (0, 2ir) interval of (f) to the range 
(0, oo) of u would mean k\ = 4, but this cannot be correct: Since the meridians 
= and = 2tt and hence u = and u = oo coincide, the corresponding rays 
must represent the same spacetime points, which is wrong, as a discussion of (|2ljj) -(|2l) j) 
shows. The correct choice is k\ = 2, which maps (0, 2tt) to the ^-interval (— oo, oo) in 
the sense that (0, tt) is mapped to (0, oo) and (tt,2tt) to (— oo,0). Since the subsets 
(v,u — > oo) and (v,u — > — oo) in the parameter space of u and v describe the same 
rays, there is no matching problem here. Taking only max regions for the v- values and 
assuming e± = 1, k 2 = \/2 for the northern and E\ = — 1, k 2 = 1/Vz for the southern 
hemisphere will satisfy our conditions. Note that the northern and southern hemisphere 
need separate copies of the max interval, u covers the range (— oo, oo), and the sign e x 
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turns out to be equal to e in Eq. (|26|). Thus the pair (9, 0) is related to a pair (u, v) by 

tan - = ( y/~( v ) \ tan — = u, (B.6) 

2 ^2 + i + ^q v 2 - v i-i) 2 v ; 

where e = 1(— 1) in the northern (southern) hemisphere. (jB.6|) may be inverted. Thus 
finally we have 

2 v^ + sinfl l-cos0 

D = — F= , u = --, 7, (B.7) 

V^-sinfl l + cos^' 

valid for the whole sphere. We can completely discard the min regions. Rays with v 
from the min region are also light cone rays, but the transformation v — > 1/v , u — > —1/u 
leads to identical rays, compare Eq. (}2T|) . thus already all rays are covered, if we 
confine the discussion to max intervals. It should nevertheless be noted that another 
parametrization is possible, which makes use of mm-intervals. 
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